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Motivated by quantum chemistry calculations, showing that molecular adsorption in graphene 
takes place on preferential sites of the honeycomb lattice, we study the effect of an isolated impurity 
on the local electronic properties of a graphene monolayer, when the impurity is located on a 
site-like, bond-like, or hollow-like position. We evaluate the local density of states (LDOS) as a 
function of energy on the impurity and on its neighboring sites, as well as in reciprocal space, at 
an energy corresponding to a bound state, in the three cases of interest. The latter study may be 
relevant to interpret the results of Fourier transformed scanning tunneling spectroscopy, as they show 
which states mostly contribute to impurity- induced variations of the LDOS. We also estimate, semi- 
analytically, the dependence of the condition for having a low-energy bound state on the impurity 
potential strength and width. Such results are then exploited to obtain the quasiparticle lifetime 
and the static conductivity in graphene in the dilute impurity limit. In particular, we recover a 
sublinear dependence of the conductivity on the carrier concentration as a generic impurity effect. 

PACS numbers: 71.23.-k, 73.23.-b, 81.05.Uw 
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I. INTRODUCTION 

Graphene is the two-dimensional allotrope of carbon, 
which is characterized by a honeycomb lattice. Despite 
its structural simplicity, only recently it has been ob- 
tained in laboratory^ 2 -! 3 -, thus giving rise to a tremendous 
outburst of research activity, both among experimental- 
ists and theoreticians. Its remarkable electronic proper- 
ties, largely due to its reduced dimensionality, and its 
relatively high degree of symmetry, make graphene an 
ideal candidate for applications in micro and nanoelec- 
tronics. In particular, it has been recently suggested that 
charging can be controlled at the atomic level, thereby 
enabling one to tailor some of the magnetic properties of 
the system^. On the other hand, its linear quasiparti- 
cle dispersion relation suggests an analogy between the 
low-energy excitations in graphene and relativistic mass- 
less particles, obeying Dirac-Weyl equation, thus allow- 
ing the study of relativistic effects in a condensed matter 
system^. 

Since most of the intriguing physical properties of 
graphene stem from its perfect crystal lattice, it is of 
interest to study how some of these are affected by the 
presence of localized impurities. It is well-known that 
disorder can significantly affect the electronic properties 
of graphene, especially when the chemical potential tra- 
verses the Dirac points. This can be brought about not 
only by impurities 7 -^, but also by topological defects^, 
edges^i, substrate corrugations^ 2 -, and ripples^ 3 -. 

Isolated short-range impurities have been shown 
to modify the local single-particle electronic proper- 
ties of graphene, such as the local density of states 



(LDOS) 7 - ^ 11 ' 14 , and can induce Friedel oscillations-. The 
role of strength, width and concentration of impurities 
in altering the local energy spectrum has been stud- 
ied theoreticall y 11 ' 15 . The relevance of special symme- 
tries and how they manifest themselves in the scattering 
around impurities has been emphasized in Ref . IrI More- 
over, the study of the impurity effects on the LDOS is 
relevant to analyze the experimental results of scanning 
tunneling microscopy (STM) 12 ' 17 i 18 , and can elucidate 
the role of correlations in the electron liquid in graphene. 
In particular, it has been suggested that Fourier trans- 
formed scanning tunneling spectroscopy (FTSTS) re- 
sults can also be instrumental to identify experimentally 
monolayer and bilayer graphene^^. 

Disorder is also known to affect considerably the trans- 
port properties of graphene. In particular, the presence 
of disorder may explain the finite value of the conductiv- 
ity in pure graphen o 20 ' 21 . In the case of graphene on a 
substrate, an inhomogeneous potential distribution may 
be brought about by charged impurities located close to 
the substrate surface. At low electron or hole concen- 
tration, this may induce sizeable spatial fluctuations of 
the carrier concentration, and may therefore justify a 
nonzero conductivity, even in the absence of any gate 
potential 21 . This has been experimentally verified using 
a scanning single-electron transistor—. Such a regime of 
inhomogeneity persists beyond neutrality, and character- 
izes also suspended graphene samples before annealing 22 . 
After annealing, the conductivity displays a sublinear de- 
pendence on carrier concentration around zero doping, 
which may be due to short-range impurities, such as point 
defects 2 ^^. 
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In this paper, we will be mainly concerned with the 
effects on the LDOS and on the conductivity of graphcnc 
due to single or distributed impurities located in a high- 
symmetry position of the honeycomb lattice. These in- 
clude the sites of the direct lattice, the position mid- 
way two neighboring carbon atoms, and the center of 
the hexagon plaquettes. Such positions have been exten- 
sively studied also within quantum chemical calculations, 
as they are expected to be favored in the adsorption of 
hydrogen, water, and other simple molecules 2 ^. 

After reviewing the formalism for a single localized im- 
purity in graphene in Sec. HU we will present our results 
for the LDOS in the presence of a single impurity, ei- 
ther in the site-like, bond-like, or hollow-like configura- 
tion (Sec. lIIIj) . Our results include the energy dependence 
of the LDOS on the impurity site and its nearest neigh- 
bors, and the reciprocal lattice structure of the LDOS 
close to a resonance. Then, in Sec. IIV1 we will gener- 
alize the above results in the case of many impurities, 
in the dilute limit, within the full Born approximation. 
In particular, we shall be interested in the case in which 
all impurities are located in a preferential class of lattice 
sites. We will derive the LDOS in reciprocal space in the 
case of many impurities, and discuss the dependence of 
the quasiparticle lifetime on the impurity concentration. 
Finally, it will be shown that, close to a low-energy res- 
onance, disorder induces a sublinear dependence of the 
conductivity on the carrier concentration, and that such 
an effect is rather insensible to the impurity concentra- 
tion, albeit in the dilute limit. We summarize our results 
in Sec. M 



II. MODEL 

We begin by reviewing the tight-binding approxi- 
mation and the T-matrix formalism for a single non- 
magnetic impurity in graphene 2 . Within the tight- 
binding approximation, a graphene monolayer in the 
presence of a single impurity localized at position x will 
be described by the Hamiltonian 



kA 



•K,* t (x)*(x). 



(1) 



Here, c kA (ckA) is a creation (annihilation) operator for a 
quasiparticle with wavevector k within the first Brillouin 
zone and band index A = 1,2, £kA = E^x — \i is the tight- 
binding dispersion relation for band A, measured with 
respect to the chemical potential fi, and Vq is a measure 
of the strength of the impurity potential. Expanding the 
field operators ^'(x), ^(x) appearing in the impurity 
potential with respect to the tight-binding basis states, 
one finds 



kA 



kk'AA' 



where 



^(k,k') = F feW4'A'(x), 



(2) 



(3) 



and ipkxfe) is the Bloch wavefunction employed in the 
tight-binding diagonalization of the pure sector of the 
Hamiltonian. 



A. Tight-binding approximation 

For the sake of completeness, we briefly review the 
main features of the tight binding approximation em- 
ployed in the present work. Graphene is characterized by 
a honeycomb lattice, with basis vectors ai = a(3, ^3)72 
and a 2 = a(3, — v/3)/2, where a = 0.142 nm is the C— 
C distance^. This is equivalent to two interpenetrating 
A and B triangular sublattices, with nearest neighbor 
sites connected by the vectors <5i = 0(1,^3)72, S2 — 
a(l,— ^3)72, S3 — a(— 1,0). Correspondingly, the first 
Brillouin zone in the reciprocal lattice is an hexagon with 
vertices in the so-called Dirac points, K = §f(L-^)j 

3a V x ' 3 /' 

A suitable choice within the standard tight-binding ap- 
proximation consists in retaining hopping and overlap 
terms between nearest neig hbor sites^. This gives rise 
to the two bands 
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(4) 



where the bottom and top signs apply to the valence 
band, with A = 1, and conduction band, with A = 2, 
respectively. In Eq. ([1J, t = 2.8 eV and s = 0.07 are 
the nearest neighbor hopping and overlap parameters, 
respectively 2 ^, and 



7k 



(•5) 



is the usual (complex) structure factor in momentum 
space. In the limit s = 0, one recovers the symmetry 
between valence and conduction bands, i?kA = ±^7k|- 

One has still a choice to fix the functional form of 
the Bloch wavefunctions that define the basis set im- 
plied in the tight binding approximation. These are lin- 
ear combination of tightly bound atomic functions, and 
will therefore be termed pseudoatomic wavefunctions in 
the following. The approximation of using pseudoatomic 
wavefunctions with a finite extension, while retaining a 
localized impurity potential, allows one to treat exactly 
also the case in which a short-range impurity is located 
in an out-of-lattice position, as is the case of hollow-like 
impurities addressed to below (Sec. IIII C[) . Due to the 
two-dimensionality of the graphene sheet, we can safely 
neglect their extension along the axis orthogonal to the 
graphene plane, z say. 

One possible choice is such that its square modulus is 
a normalized gaussiar*2£ 



(r) 



1 



2 v / 3vr a 



exp(- 



-Ps/24), 



(6) 
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where p g = Z g r/a, and Z g can be used to tune the spa- 
tial extension of the wavefunction, characterized by an 



average radius r g 



y 2 )l /2 



2V3a/Z g 



Another possible choice is such that its square modulus 
is a normalized combination of modified Besscl functions 
of second kind 31 



Mr) = -^ — ^p^ip,,) + plK ( Pb ), (7) 

where pb = Zj 3 r/a 1 and Zb is again a parameter related 
to the spatial extension. Like the gaussian pseudoatomic 
wavefunction, Eq. ©, also Eq. (JT]) has a bell-shaped be- 

3/4 

havior, but decays more slowly fa (r) ~ p b ' exp(— pb), 
for pf, ^> 1. The expectation value of any cylindrically 
symmetric function with respect to Eq. ([7]) is the same as 
the expectation value with respect to the 2p z hydrogenic 
wavefunction with atomic number Zb- In particular, the 
average radius is similarly given by f& — (x 2 + y 2 )l^ 2 = 
2\/3a/Z b . 

In both cases, the parameters are fixed by the condi- 
tion that the nearest neighbor overlap integral yields 29 
s = 0.07, so that Z g = 11.2 and Z b = 12.8. In the case of 
a single impurity, we have numerically verified that all re- 
sults do not qualitatively depend on the particular choice 
of the pseudoatomic wavefunction, and that any quan- 
titative difference is within graphical resolution. This 
is because the impurity effects considered here depend 
mainly on the short-distance behavior of (j>(r). Therefore, 
in this paper we have chosen to present results obtained 
with the gaussian choice for the pseudoatomic wavefunc- 
tions, Eq. © . In terms of these, the Bloch wavefunction 
on which the tight-binding approximation is based is 



^kA(r) 



1 



3 



(8) 



where Rj 1 are vectors of the A = A and B sublattices, 
respectively. 

Let -0k^ = A B) denote the Bloch wavefunctions 
in the sublattice representation. These are then related 
to the Bloch wavefunctions -0kA (A = 1,2) in the band 
representation, Eq. ([H]), by the unitary transformation 



fj,=A,B 

where U\^()<l) is the generic element of the matrix 



(9) 



J_( 1 -1 

V2 V~ i6k e" 



and 



7k 
l7k| 



(10) 



(11) 



B. T- matrix formalism 

We then introduce the finite-temperature Green's func- 
tions 

g A v(k,k',r) - -<T T [c kA (r)4, A ,(0)]>, (12) 

where (• • • ) is a quantum statistical average with respect 
to H at temperature T, and T T denotes ordering with 
respect to the imaginary time r. Making use of the 
fermionic Matsubara frequencies huj n = (2n + l)irkBT, 
where h is Planck's constant and fee is Boltzmann's con- 
stant, one finds the usual Dyson's equation 

t/AA'(k,k',zw„) = 5\x>Skk>G\(Kiu„) 

+ G ™' ( k > <*> *w n ) (q, k')S°< (k', iw„), (13) 

qA" 

where Q^(k, iu n ) = (iuj n ~ £kA) _1 is the Green's function 
of the pure system. Dyson's equation (|13[) can be solved 
iteratively by exploiting the fact that the impurity poten- 
tial is factorizable in momentum space (see Appendix [A"]l . 
One finds 

£?AA' (k, k', ito n ) = 6\y6wGl(k,iuj n ) 

+g° x (k, iu n )T xx , (x; k, k', kj n )g° x , (k', iu n ), (14) 



where 



T A a'(x; k, k.',iuj n ) 



1 ^yxjj^x) 
7V1-V £°(x,x, iw n ) 



(15) 



is the generic element of the T- matrix, ^^(x) = 
V^V" ^kA ( x ) is a rescaled basis function, and 

g°(r,r',w„) ^^^tw^M. (16) 

qA' 

Eq. (fTi|) shows that the correction due to a single lo- 
calized impurity vanishes as l/N in the thermodynamic 
limit, N — > oo. Going back to real space by means of 
Eq. (|16p . one finds the imaginary-time Green's function 
at position r 



Vbg°(r, x, zu n )g°(x, r, joy, 
1 - V G°(x,x, iu n ) 



g(r,r,iuj n ) = g°(r,r,iu} v 

! - I „y (A. A. I 

(17) 

In what follows, we shall be interested in the local density 
of states (LDOS) p(r, oj), which is experimentally accessi- 
ble through STM measurement o 12 ! 17 ' 18 , and is related to 
the imaginary part of the analytically continued Green's 
function through 



p(T,w) = --ImG(r,r,w), 



(18) 



with 7k defined in Eq. (0. 



where G(r, r, oj) = g(r, r, ioj n — > ui + i0 + ). It is straight- 
forward to observe that the chemical potential enters 
G(u>) only as an additive constant to u>. Therefore, we 
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(a) site-like 
(b) bond-like, x 5 
(c) hollow-like, x 2 10 



rocal space, which reads 




(a) site-like 
(b) bond-like, x 5 
(c) hollow-like, x 2 10 




FIG. 1: (Color online) Real part of the unperturbed Green's 
function, ReG°(x,cj) (top panel), and unperturbed LDOS, 
p°(x, lo) = — 7r _1 Im G°(x,a;) (bottom panel), for the three 
cases of interest: (a) site-like impurity (x = 0), (b) bond-like 
impurity (x = S3/2), (c) hollow-like impurity (x = S3). It 
should be noticed that the latter two cases have been scaled 
by the factors indicated in the caption. 



can set hereafter p = 0, thereby neglecting any contri- 
bution arising from chemical or electrical doping, e.g. 
through a gate voltage. 

Inspection of Eqs. $F7§ and (18|) shows that the LDOS 
on the impurity position (r = x) is given by 



p°(x,c) 



[1 - VoRe G°(x, x, w)]a + [ttVqp ^ u)f 



(19) 

with p°(x,lo) denoting the LDOS at position r = x 
in the pure case. In the limit of a vanishing unper- 
turbed LDOS, p°(x, to) -> 0, one has p(x, to) -> V ~ l S[l - 
VoReG (x, x, ui)]. Such a circumstance is e.g. realized 
below the valence band, lo < lo± w — 2.48i, or above the 
conduction band, lo > ojj ~ 3.80t. Direct inspection of 
Re G°(x, x, lo) as a function of lo (Fig. [T]) leads to the ex- 
istence of a bound state outside the two bands, for a wide 
range of potential strengths Vq. In particular, a bound 
state at lo < may be formed below the valence band 
only for some Vq < 0, or above the conduction band for 
Vq > 0, the energy of such a bound state moving farther 
from the bands, as |Vb| increases. 

For future reference, it is also of interest to quote the 
expression of the LDOS close to a bound state in recip- 



PA(k, w) 



1 



— ImGA(k, lo) 



Vo \M*)\ 2 
N |^-£ Ak | 2 



6[l - VbReG°(x,x,w)](20) 



In other words, Eq. ([20]) applies to states with a vanish- 
ing unperturbed LDOS, i.e. to frequencies lo such that 
p°(k, lo) = for all wavevectors k in the 1BZ. This cor- 
responds to u> < lo±, lo — 0, and lo > lot- 



III. SINGLE IMPURITY 

In what follows, we shall analyze the effect on the 
LDOS, Eq. (fl8)) . due to a single impurity localized in 
several high-symmetry positions of the primitive cell in 
the graphene honeycomb lattice. These include an A or 
B site, usually occupied by a carbon atom (site-like im- 
purity), the position midway between an A and B site 
(bond-like impurity), and the position at the center of 
an hexagon plaquette (hollow- like impurity). 



A. Site-like impurities 

Let us start by considering an impurity located on 
an A or B site (say x = 0, for definiteness) . Such an 
impurity preserves the D^h symmetry. This could be 
used to model a hydrogen impurity adsorbed on a carbon 
ato m 16 ' 26 , or a vacancy (here obtained in the Vq — > 00 
limit) 26 ' 32 , as could be induced by proton irradiation 3 ^. 

Fig. [2] shows the LDOS on a single site-like impurity, 
Eq. (Hn|), for negative as well as for positive values of 
the impurity strength Uq — Vq/o 2 , where a is the lat- 
tice step, in the limit of weak scattering (|?7o| 6i). 
As anticipated, a bound state forms below the valence 
band if Uq < 0, whereas a bound state forms above the 
conduction band for Uq > 0. 

From Eq. (|19[) . a resonance is formed at an energy Lo res 
between the two Van Hove singularities when 



1 - VbReG (x,x,w res ) = 0. 



(21) 



By inspection of the w-dependence of Re G°(lo) (Fig. [T]), 
it follows that Eq. (JH]) is fulfilled for -1 < t/U < \. 
Such a resonance is better resolved when the unper- 
turbed LDOS /0°(x, x, lo) is weak for lo w ui res . This 
is indeed the case in the proximity of the Dirac points, 
where p°(x,x, lo) — > linearly as lo — > 0. A special case 
is represented by the limit Uq — > 00, corresponding to a 
vacancy formation. In this case, the condition for a res- 
onance is fulfilled at lo m st, where however the LDOS is 
strongly depressed. 

At exactly u> = 0, e.g. when the chemical potential 
traverses the Dirac points, the resonance becomes a true 
bound state, since /?°(x, x, to) = 0. The value of the im- 
purity potential Uq allowing a bound state at lo = can 
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FIG. 2: (Color online) Local density of states p(x = 0, u>), 
Eq. (llQfl . on a site-like impurity located at x = 0. Top panel 
shows the LDOS for U /t = -0.05, -0.10, -0.15, -0.20. Bot- 
tom panel shows the LDOS for U /t = 0.05, 0.10, 0.15, 0.20. 
In both panels, we also show the LDOS in the unperturbed 
case (Uo/t = 0). 
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FIG. 4: (Color online) Showing the LDOS for a site-like impu- 
rity at x = (a) on the impurity site, (b) on a nearest neigh- 
bor site, (c) on a generic lattice site in the unperturbed case. 
Top panel refers to low potential strength (Uo = Lit), while 
bottom panel refers to large potential strength (Uo = 3.7t). 




FIG. 3: (Color online) Contour plots of the LDOS in momen- 
tum space, p\(k, lo = 0), Eq. (|20p . for the valence (A = 1, 
left panel) and conduction band (A = 2, right panel). Here, 
we are considering a site- like impurity with Uo = 3.7i, thus 
giving rise to a bound state at lo — 0. 



be obtained within a semi-analytical approach (see Ap- 
pendix [B] for details) . This is based on an expansion of 
ReG°(x,x,u = 0) in Eq. flU} at x = 0. To the lead- 
ing terms in the nearest neighbors, one finds that the 
condition for having a bound state at lo = is 

±tt<j>(0)[scf>(0)+2 ( f>(S 1 )}, (22) 
Uo 



where </>(x) is a gaussian pseudo-atomic wavefunction, s 
is the band asymmetry parameter, and the position 
of a nearest neighbor to x = (Appendix [B| . The first 
contribution to Eq. (122p is due to the band asymmetry 
(s 7^ 0), whereas the second contribution is related to 
the wavefunction width, and can be neglected for a suf- 
ficiently localized pseudo-atomic wavefunction. In the 
limit of symmetric bands (s = 0) and localized wave- 
functions, one recovers a bound state at exactly lo = 
in the case of a vacancy (Uo = oo) ^ 11 ' 15 . From Eq. |22|) 
one may conclude that a bound state is formed also in the 
case of a localized pseudoatomic wavefunction, provided 
one retains a nonzero band asymmetry (s ^ 0), and that 
this takes place for a finite value of the impurity poten- 
tial (Uo < oo ), in agreement with the findings of Ref. 7. 
Fig. [3] shows a contour plot of the LDOS in momentum 
space for an impurity potential generating a bound state 
at lo = 0, for both the valence and conduction bands. 
In both cases, the largest contribution to p\(lt, lS) comes 
from the wavevectors close to the Dirac points. Slight 
differences between the two bands are due to a nonzero 
asymmetry parameter s. 

We end this subsection by considering the effect of a 
site-like impurity located at x = on the LDOS at a 
neighboring site, y = S3, say. After the appropriate ana- 
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lytical continuation, Eq. (JTTj) then yields 

/ s ReG (y,x,u,)ReG°(x, y ,^) 

Pyy^l ~ rp n0 , via 

[ReG u (x, x, io)\ z 

xF " 1 5[l-y ReG (x,x,w)], (23) 

in the limit of vanishing unperturbed LDOS on the im- 
purity, p°(x, ui) — > 0. Therefore, while the condition for 
having a bound state on a neighboring site is the same 
as Eq. (|2"Tj) , the relative weight with respect to the im- 
purity site is given by the prefactor in Eq. (|23|) . Fig. [4] 
shows the LDOS on a neighboring site (y = 63), when 
a site-like impurity is located at x = 0. The LDOS cor- 
responding to a bound state outside the band is larger 
on the impurity site than on the neighboring site. The 
opposite is true for the LDOS corresponding to the reso- 
nant state between the two Van Hove singularities, which 
is depressed on the impurity site than on the neighboring 
site. The same effect applies to the bound state between 
the two Van Hove singularities. 

Such a behavior for resonant states in the energy range 
between the two Van Hove singularities is analogous to 
the one encountered in the d-density-wave (DDW) phase, 
which has been suggested as a viable description of the 
pseudogap phase in the high-T c cuprate a 34 ' 35 . In the 
DDW phase, the LDOS on an atomic site vanishes lin- 
early as 00 — ► for the pure system, and exhibits two 
Van Hove singularities, symmetric with respect to fi = 0. 
It has been demonstrate d 36 ! 37 that a sufficiently strong 
localized impurity produces a resonance between the two 
Van Hove singularities. The LDOS associated to such 
a resonance in the DDW phase is larger on a nearest 
neighbor, than on the impurity site, in close analogy to 
what is here shown for an impurity in graphene, and in 
agreement with the findings of Ref. [Ij. In both cases, the 
quasiparticle bands are characterized by two inequivalcnt 
minima (the two Dirac points, in the case of graphene), 
so that scattering processes due to short-range impurities 
can be classified as intra-valley or inter-valley, depending 
on whether initial and final states lie close to the same or 
to different extrema, respectively. The peculiar behavior 
of the LDOS corresponding to resonant states at u> « is 
related to inter-valley scattering, at variance with bound 
states outside the bandwidth. It is also relevant to note, 
in this context, that single-impurity scattering around a 
localized impurity has been suggested as a tool to distin- 
guish between a DDW and a pseudogap phases, within 
the precursor pairing scenario 38 . 

B. Bond-like impurities 

An impurity located between an A and B site only 
preserves the C<2. v symmetry, and may be used to model 
an oxygen impurity between two carbon atoms^ 6 -. This 
corresponds to three inequivalent positions in the real 
lattice, although local effects on each of them are related 
by rotations of multiples of 27r/3. In the following, for 




-3-2-1012 3 4 5 -3 -2-1 1 2 3 4 5 



FIG. 5: (Color online) Local density of states for a site-like 
impurity (top row panels), and a bond-like impurity (bottom 
row panels). Left panel on top refers to the LDOS evaluated 
on the same position of a site- like impurity (r = x = 0); 
right panel on top refers to the LDOS for a site-like impurity 
(x = 0), evaluated on a nearest neighbor lattice site (r = Si). 
Left panel on bottom refers to the LDOS evaluated on the 
same position of a bond-like impurity (r = x = 53/2); right 
panel on bottom refers to the LDOS for a bond-like impurity 
(x = 83/2), evaluated on a nearest neighbor lattice site (r = 
0). The potential strengths are (a) Uo = 0; (b) Uo = 0.05(7o; 
(c) Uo = 0.10f7o; (d) Uo = 0.25(7 ; (e) U = 0.50£/ ; (f) 
Uo = 0.75U ; (g) Uo = U ; (h) U = 2.00(7 , where U is the 
value of Uo yielding a bound state at uj = 0. 

definiteness, we shall therefore be concerned with a bond- 
like impurity located at x = 63 /2. 

With reference again to Eq. (fT9|) , one finds a markedly 
different w-dependence of G°(x, x, to) at x — 83/ 2 
(Fig. [TJ, with respect to the site- like case (x = 0). In- 
deed, while ImG°(w) is depressed with respect to the 
site-like case, as a consequence of the finite extent of the 
gaussian pseudoatomic wavefunctions, ReG°(w) attains 
a finite value at uo = 10^ (Fig. [1]). This implies the ex- 
istence of bound states above the conduction band only 
for Uo > 8t, while bound states below the valence band 
still exists for any Uo < 0. On the other hand, resonances 
between the two Van Hove singularites close to 00 = are 
possible for 2 < Uo/t < 20, i.e. within a finite range of 
positive values of the impurity strength. 

An expansion of ReG° in Eq. (f2Tj) . where now x = 
S3/2, leads to the estimate 

±^(A b + 2 s )cb 2 (S 3 /2) (24) 
Uo 

for the impurity potential required to generate a bound 
state at 00 = 0, where ~ 0.67 (Appendix [B"|) . As in 
the site-like case, one can recognize a term due to the 
asymmetry between the two bands (s ^ 0). 

Fig. O compares the LDOS for a site-like (top row 
panels) and a bond-like (bottom row panels) impurity, 
both evaluated on the same position as the impurity 
(left column panels) and on a nearest neighbor lattice 
site (right column panels), for several values of the po- 



FIG. 6: (Color online) Contour plots of the LDOS in mo- 
mentum space, p\(k,u>), Eq. (|20[1 . for the valence (A = 1, 
left panel) and conduction band (A = 2, right panel). Here, 
we are considering a bond-like impurity with Uo ~ 5. It, thus 
giving rise to a bound state at lo — 0. 



tential strength. One finds quite a different behavior in 
the two cases. As Uo increases towards Uo, i.e. the 
value of the impurity strength giving rise to a bound 
state at to = 0, the LDOS on the impurity site, and the 
weight of the Van Hove singularities, decreases, as ex- 
pected. On the other hand, the bound state at w = 
becomes sharper and more pronounced in the bond-like 
case, whereas it becomes suppressed in the site-like case. 
Another remarkable difference is the presence, in the 
bond-like case, of a wide resonant state in the conduction 
band for 4t < Uq < 8i, which is completely absent in the 
site-like case. This can be traced back to the different in- 
dependence of Re G° in Eq. (|2ip in the two cases. Fig. [5] 
also compares the LDOS for the site-like and bond-like 
cases, but now evaluated on a nearest neighbor lattice 
site to the impurity position. Again, the Van Hove sin- 
gularities become smoother, as Uq increases towards Uo. 
The suppression of the singularities due to an increase in 
the potential strength is enhanced in the bond-like case, 
than in the site-like case. 

Finally, Fig. [S] shows the LDOS in momentum space, 
Eq. (pQ|). in the case of a bond-like impurity at x = S 3 /2, 
for the valence and conduction bands. Fig. [6] refers to 
a potential strength of Uq = 5.1t, thus giving rise to a 
bound state at ui = 0. Similar pictures, but rotated of 
multiples of 27r/3, would be obtained in the other, in- 
equivalent, bond-like positions. As in the site-like case, 
Fig. [3l one finds that the points in k-space providing 
the largest contribution to p\ (k, ui) are those closer to 
the Dirac points, but now with a reduced symmetry. In 
particular, p\(\t, ui) is not invariant with respect to trans- 
formations of the Cq v point group because of the squared 
modulus of the pseudoatomic wavefunctions in Eq. ([2U)l . 



C. Hollow-like impurities 

The last case considered here corresponds to having 
a single impurity located at the center of an hexagon 
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FIG. 7: (Color online) Local density of states for a hollow- 
like impurity (x = —S3) on the same site as the impurity (top 
panel, r = x), and on a nearest neighbor lattice site (bottom 
panel, r = 0). The potential strengths are (a) Uo = 0; (b) 
Uo = 0.10f7o; (c) Uo = 0.25C/ ; (d) Uo = 0.50C/o; (e) Uo = 
0.75Ub;_(f) U = Uo; (g) U = 1.50U ; (h) U = 2.00t/ , 
where Uo is the value of Uo yielding a bound state at lo — 0. 

plaquette, x = — S3, say. This is the highest symmetry 
position in the carbon honeycomb lattice, and indeed the 
point symmetry D 6 ?i is preserved. Inspection of Fig. ^ for 
the w-dependence of ReG°(w) and p°(ui) shows that the 
unperturbed LDOS is severely depressed (some three or- 
ders of magnitude lower) than the LDOS in the site-like 
case. Moreover, as a consequence of the overall behav- 
ior of ReG°(w), one has a bound state below the va- 
lence band for all negative values of Uo, whereas one has 
a bound state above the conduction band for relatively 
large positive values of the impurity strength, Uo > 10 4 t. 
On the other hand, the relatively low value of p°(uj) al- 
lows the formation of well-resolved resonant states close 
touj = 0, for 1.5-10 3 i < U < 1.1-10 4 *. Expanding Re G° 
in Eq. (|2ip. where now x = — 83, yields in this case the 
estimate 

4- « (A h + 2B h s)<t> 2 (6 3 ) (25) 
Uo 

for the impurity potential required to generate a bound 
state at ui = 0, where Ah — 2.35 and Bh = 3 (Ap- 
pendix [5]). As in the previous cases, the main term per- 
sists also in the limit of perfect band symmetry (s = 0). 

Fig.[7]shows the LDOS for a hollow-like impurity, both 
on top of the impurity site, and on an adjacent lattice 



8 




FIG. 8: (Color online) Contour plots of the LDOS in momen- 
tum space, PA(k, lu), Eq. (|20|l . for the valence (A = 1, left 
panel) and conduction band (A = 2, right panel). Here, we 
are considering a hollow-like impurity with Uo = 3.8 ■ 10 4 t, 
thus giving rise to a bound state at lu = 0. 



site, for several potential strengths. At variance from 
the previous two cases, it is apparent that resonant states 
are sharper in the nearest neighbor site, than on top of 
the impurity position. Analogously to the bond-like case, 
there are resonant states developing in the high conduc- 
tion band, which are however better defined. 

Fig. [8] shows the LDOS in momentum space, for a 
hollow-like impurity giving rise to a bound state at u> = 0, 
Eq. (|25p . At variance with the previous cases, one may 
notice that the k-states contributing most importantly 
to p\(\i, lu) are the same as those involved in building up 
the unperturbed LDOS. Indeed, in the conduction band 
(u) > 0), the largest contributions to p°(— S3, — S3, lu) 
come from the Van Hove singularities and the centers 
of the sides of the first Brillouin zone. Similarly, in 
the valence band (w < 0), the largest contributions to 
p°(— S3, — S3, lu) come from the band bottom, i.e. from 
k-points close to the T point. This can be ultimately be 
traced back to the extended width of the gaussian pseu- 
doatomic wavefunction here employed. 



IV. MANY IMPURITIES 



fects can be neglected; (iii) their number is sufficiently 
large (N imp ^> 1), so that their effect is appreciable on 
bulk properties in the thermodynamic limit, but the im- 
purities are sufficiently diluted (ni mp = Ni mp /N -C 1). 
In these limits, while a standard averaging procedure over 
the position configurations of the impurities restores the 
translational invariance of the Green's function, 



^7(k,k',^ n ) = 5 kk ^i7(k,io;„), 



(26) 



the eigenstates of the pure Hamiltonian are expected to 
acquire a finite lifetime induced by disorder. This can be 
formally achieved by relating Q lmp (k, ioj n ), now a matrix 
with respect to the band indices, to the single- impurity 
Green's function Q°(k,iuj n ) discussed in Sec. IIIII through 
a Dyson's equation analogous to Eq. (fT3|) , but now in- 
volving the proper self energy matrix^ 7 . S(k, iu> n ) 

g m p(k,iu; n ) = g°(k,iu n ) 

+G°(k, io; n )E(k, iLu n )G tmp (k, iLu n ). (27) 

Within the full Born approximation (FBA)2£, which is 
valid in the limit of small impurity concentration, n? mp -C 



-XX 



(k, iu> r , 



Vb V 
N i-i 



q,\ 



» |V'kA"(x)| 2 ^„(q,iw„) 

where we are assuming that all impurities occupy equiva- 
lent lattice positions. Comparing Eq. (|28f with Eqs. (fT5|) 
and (fT6| . it is possible to relate the proper self-energy 
within the FBA to the T-matrix for the same kind of 
impurity, through 



J AA' 



N im pT\\>(Tf,'k,'k',iuj r , 



(29) 



Eq. (|29|) therefore enables us to generalize most of the 
results derived in Sec. IIIII to the case of many impurities, 
all located within a preferential class of lattice positions. 



While single impurity effects are in principle observ- 
able through STM measurement a 12 ' 17 ! 18 , real samples 
usually contain a sizeable amount of impurities, which 
are responsible of sensible modifications of both thermo- 
dynamic and transport properties. Therefore, we will 
here exploit the results of Sec. IIIII for a single impurity, 
to study the effect of Ni mp impurities on a graphene 
monolayer. We will assume that (i) the position of all 
impurities differ by a vector of the direct lattice; in 
other words, there is a preferential kind of impurity lo- 
cation, i.e. all impurities are either site- like, bond- like, 
or hollow-like, according to the classification given in the 
Sec. IIIII (ii) impurities are independent, i.e. the aver- 
age distance between two impurities is larger than the 
quasiparticle coherence length, so that interference ef- 



A. LDOS 

We begin by discussing the effect of many impurities on 
the local density of states in reciprocal space. This can 
be obtained by inverting Eq. I|27p for g trnp {k, iu> n ) and 
then performing the usual analytical continuation. Most 
properties can be derived by describing the behavior of 
the analytically-continued proper self-energy, which in all 
of the three cases of interest can be written as 



£(k,w) = 



V W(k) 



L imp 



l-Vo<3°(x,x,u;) 



(30) 



where x = 0, S3/2, or — S3 in the site-like, bond-like, or 
hollow-like case, respectively, and VF(k) is a matrix form 
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factor explicitly given by 



W$(k) = ifc(0)4v(0)+^ kA (^)^kv(<53ffla) 



1 3 

~E^a(<W2)^v(<V2) (31b) 



^kA(-<53)^kA'(-<53) (31C) 



in the site-like, bond-like, and hollow-like cases, respec- 
tively. We are here assuming that the impurities are 
equally distributed among the A and B sites, in the site- 
like case, and among the three classes of a bonds, in the 
bond-like case. 

Both in the site-like and in the bond-like cases, di- 
rect inspection of the solution of Eq. (f2"T)) shows that 
G lmp (k, oj) is nearly diagonal in the diluted limit (n 2 mp <C 
n-imp)- Therefore, an eigenstate of the unperturbed 
Hamiltonian labelled by wavevector k and band index 
A acquires a finite lifetime t^a, which e.g. in the site-like 
case and in the limit of low LDOS is given by 



'kA 



'imp 



V 2 W xx (k)p(0,Lo = & x ), 



(32) 



where p(x, lu) is the LDOS with a single impurity, 
Eq. (HD. 

For impurity potentials close to the condition for a 
well-defined resonance at lu ss in the single-impurity 
case, Eq. (|2"TT) . in the dilute limit, one gets for the LDOS 
in reciprocal space close to a Dirac point 



PA(k,w) 



.4 



kA 



k-£kA + £kA(^-£kA)] 2 + ^kA' 



(33) 



where A = 2 when A = 1, and A = 1 when A = 2, and 



case, one finds a nearly diagonal Green's function, whose 
nonzero matrix elements are given by 



Gimp 
AA 



W xx (k) f 
W-aaW 

{ i W xx Q<.) (u-t Xk ) 2 
7T Wf x (k) n imp V^p(-5 3 , uj) 



(35) 



In the same limit, for an impurity potential giving rise to 
a resonance at exactly lu = 0, the situation is even more 
dramatic, since 



°AA 



-u 



;Ak 



V^AA(k) 

W xx (k) 



Ak 



(36) 



where r\ is a positive infinitesimal. On the other hand, 
#3, CiVes)|, one finds 



for n imp <§; t |ImG°( 



Gimp 
XX 



- 6>k + i7rn imp V 2 p(-(53,a;)W / AA(k)] 



(37) 

The different behavior with respect to the previous two 
cases follows from the fact that in the hollow-like case the 
Born approximation performs an average with respect to 
impurity positions all of the same kind, at variance, e.g. 
with the site-like case, where impurities can be added 
either in the A or in the B sublattices. This is likely 
to produce an additional dephasing among contributions 
arising from inequivalent lattice positions in the site-like 
or bond-like cases, with respect to the hollow-like case, 
thereby resulting in an increasing inverse lifetime with 
increasing impurity concentration, although only in the 
moderately low impurity concentration, Eq. (|35|) . One 
however recovers the physically expected behavior at low 
impurity concentrations, i.e. a vanishing inverse lifetime 
with increasing impurity concentration, Eq. (|37p . 



.4 



kA 



-BkA — 



n imp V<?W\x(k)px(0,u) 
Wi 2 (k)W 2 i(k) 

Wn(k)W22(k) 
W 12 (k)W 21 (k) 



(34a) 
(34b) 



The behavior of the LDOS in reciprocal space is therefore 
quite different from the unperturbed case, which would 
be characterized by a Dirac delta peaked along closed 
contours around the Dirac points. 

The case in which all impurities are located in a hollow- 
like position, Eq. (|31c[) . is quite different from the previ- 
ous two cases, Eqs. (I31al) and (|31b[) . This is due to the 
fact that the form factors VK(k) defined in Eq. (|34c|) fulfill 
the additional identity Wii(k)W22(k)-Wi 2 (k)W2i(k) = 
0. 

While the Born approximation holds for low impurity 
concentration, n 2 mp -C ni mp , one has to distinguish two 
regimes. For moderately large impurity concentrations, 
nimp S> t |ImG °(— #3,u; res )|, for an impurity potential 
close to what would be a resonance in the single- impurity 



B. Conductivity 

We end this section by considering the effect of many, 
short-range impurities on the conductivity. As described 
above, we are mainly concerned with the case in which 
all impurities are located in the same class of lattice po- 
sitions, and we will here focus on the site-like case, even 
though, as described in Sec. IIV Al the analysis presented 
here is actually more general, as it qualitatively applies 
also to the bond-like case, at least for u between the two 
Van Hove singularities. Again, we will assume the dilute 
regime, n 2 mp <C Ui mp , so that the full Born approxima- 
tion holds. 

Within linear response theory, the conductivity a is re- 
lated to the current-current correlation function through 
a Kubo formula 



<7i m (u,T;u>) 



-Sir; 



hujNA ce ii 



IlL (0,0,u,), (38) 



where n is the electron density, u> is the frequency of the 
external electric field, A ce u is the area of a primitive cell, 
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and nJJ n (k, q, us) is the (l,m) component of the Fourier 
transform of the retarded current-current correlation ten- 
sor. We are mainly interested in the dissipative part of 
the conductivity tensor, i.e. its real part. For the longi- 
tudinal part a — Re a xx , one has 

a(^) = -^J^lmflt(0M, (39) 

where nj^ is the retarded version of 

fl /m (k,q,r) = -<T T [j^(k, r)j^ (q,0)]>, (40) 



and Jj V (k, r) denote the Fourier transform of the param- 
agnetic component of the current density vector, at the 
imaginary time r. Expanding J ; v (k, r) in terms of the 
eigenstates of the unperturbed Hamiltonian, and making 
use of the results of Appendix [Cl one finds 



11^(0,0^)= e 2 -^- ]T ^, A ,(k)/ lx ,v,'(k')(T r [ C l cA (r+) Ck ,(r)4, A ,(0 + ) Ck , t ,K0)]>, (41) 

kk' 

where the matrix elements h m .aij(x) are defined in Appendix [Cl One is now in the position to make use of Wick's 
theorem. We further make the approximation to treat the one-body Green's functions within the FBA, and perform 
the required analytical continuation, to obtain 

y-2 2 poo ,f 

n*(o,o,") = ^— VA*(i0^v(k)y_ M ^{M^ 

+MC)G;T(k, OG^k, C + w) - n F (C + ^)G A 7*(k, C)G™f(k, C + w)} , (42) 

I 



where np(w) is the Fermi function, and the factor of 2 
takes into account for spin degeneracy. 

We next make a further approximation, i.e. we assume 
that the impurity Green's functions are diagonal in the 



^AA'G™ p (k, w). Such an ap- 



band index, G A ™f (k, a;) 
proximation is justified in the dilute limit, and amounts 
to treat the effect of disorder as a perturbation to the 
pure case, whose main effect is that of adding a finite 
lifetime to the eigenstates of the unperturbed Hamilto- 



In the static limit (u> — > 0) and at T = 0, Eq. 
yields the conductivity as a function of the chemical po- 
tential, which can be decomposed in an inter and intra- 
band contribution, <tdc(a*) — Winter (m) + (fi), given 

by 

amtM = -^£^,i2(k)Wk) 



x Im G l ™ p (k, 0) Im G™ p (k, 0) (43a) 



44 E [^, QQ (k)ImG»S P (k,0) 



CTQ 



k,a = 1.2 



(43b) 



where Co = 7re 2 /(2/i) is proportional to the quantum 
of conductivity, and t q ~ 2 = 16t 2 /(3V3nh 2 ). One may 
expect that the interband contribution, Eq. I|43a|) , only 



becomes comparable with the intraband contribution, 
Eq. (|43b[) , when /x sa 0, i.e. when the valence and conduc- 
tion bands overlap, owing to the disorder-induced energy 
spread. Away from neutrality (/x = 0), and for a given 
impurity potential strength Uq, an increase of the Fermi 
surface width produces an increase of the conductivity. 
Such an increase is however rather slow, close to the en- 
ergy values where the LDOS with a single impurity is 
maximum, where the backscattering due to the impuri- 
ties is more effective. Such a sublinear increase of ctdc as 
a function of the carrier concentration n occurs for values 
of Uq giving rise to resonant states close to fi = and 
does not depend on the value of rii mp . Such a behavior 
is numerically confirmed in Fig. [9l for various values of 
n%mp and Uq, and is in good qualitative agreement with 
the experimental results^ 2 -. The asymmetry between the 
particle (n > 0) and hole (n < 0) regimes is partly due 
to the band asymmetry (s ^ 0), but is mainly due to 
the effect of impurities, which is different depending on 
the sign of //. Both in the valence and conduction bands, 
however, one observes the occurrence of a maximum and 
then a decrease of odc when /i attains the value cor- 
responding to a Van Hove singularity, where the Fermi 
surface is maximally extended and traverses an electronic 
topological transition. A comparison between the two 
panels in Fig. [9] shows that the nonmonotonic behav- 
ior of ctdc is generic for all impurity concentrations, but 
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rather depends on the impurity potential Uq. The simi- 
larity between Fig. [5] and the concentration dependence 
of the conductivity measured in suspended graphene af- 
ter annealing surmises that scattering due to short range 
impurities is relevant to determine the transport proper- 
ties of these graphene sample o 23 ' 24 ^ 39 . 



V. CONCLUSIONS 

We have analyzed the effects of a single, localized im- 
purity on the local electronic properties of a graphene 
monolayer. Specifically, we have considered an iso- 
lated impurity located on high-symmetry positions of 
the honeycomb lattice, such as the site-like, bond-like, 
and hollow- like positions. While the electronic proper- 
ties of the pure system have been treated within the 
tight-binding approximation, but allowing for asymme- 
try between valence and conduction bands, the effect of 
the impurity has been modelled through a gaussian pseu- 
doatomic wavefunction, even though more general func- 
tional forms have been taken into account. Moreover, the 
tight binding scheme employed in this work does not suf- 
fer from the 'cone approximation ' 11 ! 15 , thereby enabling 




one to treat both low and high energies within the band 
with the same degree of accuracy. 

We have evaluated the local density of states as a func- 
tion of energy on the impurity site and on its nearest 
neighbor locations, and as a function of wavevector in 
reciprocal space, close to a resonance. The latter may be 
of relevance to interpret FTSTS measurements around an 
impurity. In particular, it has been shown that the main 
contributions to the impurity- induced modification of the 
LDOS come from wavevectors close to the Dirac point, 
in the site-like and bond-like cases, while the same states 
are involved in determining the LDOS both in the un- 
perturbed case and in the case of a hollow-like impurity. 
Moreover, it has been suggested that FTSTS spectra can 
be used to distinguish between different types of impuri- 
ties, in particular as far as the impurity potential range 
is concerned^. 

We have determined semi-analytically the condition on 
the impurity strength for having a bound state at to = 
in the three cases of interest. In particular, in the site-like 
case, it is shown that the weight associated with a bound 
state between the two Van Hove singularities is larger for 
the LDOS on a neighboring site, than on the impurity site 
itself. Such a behavior is analogous to the one predicted 
for the d- density- wave state of high-T c superconductors, 
and can be attributed to the different contributions com- 
ing from intra- and inter-valley impurity scattering. 

Our results for the single-impurity case have been ex- 
ploited to discuss the effect of distributed impurities, all 
located in a preferential class of lattice sites. Such a gen- 
eralization has been derived within the full Born approx- 
imation, and applies to the dilute limit. In particular, 
we have estimated the quasiparticle lifetime associated 
to a finite impurity concentration, and the behavior of 
the LDOS in reciprocal space. Within linear response 
theory, we have also evaluated the static conductivity. 
One can again distinguish an intra- and an inter-band 
contribution, the latter being sizeable only close to zero 
carrier concentration, i.e. when the two bands apprecia- 
bly overlap. Moving away from neutrality, one recovers a 
nonmonotonic dependence on the carrier concentration, 
characterized by a sublinear increase close to fi — 0, as is 
observed experimentally in suspended graphene samples 
after annealing. Such a feature is generic, in the sense 
that it applies to all impurity concentrations under study, 
and rather depends on the impurity strength. 
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FIG. 9: (Color online) Conductivity <jdc as a function of the 
carrier concentration n in graphene, in the presence of many 
site-like impurities. The upper panel refers to a value n; mp = 
1CP 2 of the impurity concentration, while the lower panel is 
characterized by ni mv = 10~ 3 . The impurity potentials under 
consideration are (a) Uo = 0.35?7o, (b) Uo = 0.50£/o, (c) Uo = 
Uo, (d) Uo = 10.00(7o, (e) Uo = — Uo, where Uo is the potential 
strength giving rise to a bound state at to = 0. 
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APPENDIX A: DYSON EQUATION FOR 
SEPARABLE IMPURITY POTENTIAL 



Here, we briefly derive Dyson's equation, Eq. (|T4|) . 
in the case of a separable impurity potential, Eq. ([3]). 
Inserting our Ansatz for the potential, Eq. ([3]), into 
Eq. (| 13(1 . and iterating, one may express the Green's func- 



tion Q\\i (k, k', icu n ) as a series*! 



e?Av(k,k',io;„) = ^ail ) ,(k,k / ,ia; n ), (Al) 

1=0 



whose first and successive terms are given iteratively by 



gW(k,k',iu n ) = 5 X y5 kk ,g { ^(k,iLu n ) (A2a) 
G[% (k, k', iw n ) = gfx l) (k, k', iun) + V J2 Sxy^ (k, q, WVk'A' (x)^ (k', iWn). (A2b) 



7 (0), 



qA" 



This leads to the series 

gw (k, k', iuJn) = <5AA"5kk'^i 0) (k, iw„) 
?(°)/ 



+K,^ u ^(k, z Wn )^ A (x)Vk'v(x)^? ) (k', iw„) 



+K ) ^ 0) (k,zc n )V k A(x) ( 2Vo<A«W??(q,»Wn)V'qA«(x) I VVA'WSi^k',^™) 

, qA" 



+K ) £ A ° ) (k, W„)fe(x) ( ^ VbV'qA" (q, iw n )VqA« (x) ] fax- (x)^? } (k', M „) 

> qA" 



which is recognized as a geometric scries, whose sum can be cast in the form of Eq. (114 

| 



(A3) 



APPENDIX B: EXPANSION OF THE BLOCH 
WAVEFUNCTIONS IN THE SUBLATTICE 
REPRESENTATION 

At the origin of the A sublattice, r = say, the Bloch 
wavefunctions can be expanded as 



^kA(O) 
^kB(O) 



J A ^k 



J A ^k 

t(2) (2) 
^ 7k 



(Bla) 
(Bib) 



where (3^ and 7^ are basis functions of the trivial irre- 
ducible representation of the point group D§h and -D3/1, 
respectively. In particular, one finds = 7k, Eq. (0, 
while = /3k, with 



'k 



(B2) 



where the prime restricts the summation to all next near- 
est neighbors in the direct lattice, i.e. |Rj| = V^a. 

Because of the rapid decrease of the gaussian pseu- 
doatomic wavefunction, Eq. (JSJ), one may safely truncate 
the expansions Eqs. (|B1[) to the first terms, thereby ob- 



taining 



^kA(O) 

V>ks(0) 



N 



•K^ihk 



(B3a) 
(B3b) 



The latter can be used in the expansion of 
ReG°(x,x,w = 0) appearing in the resonance con- 
dition, Eq. (|2ip . which for a site- like impurity reads 



ReG°(x,x,0) = -]>>kA(x)| 2 

kA 

+Vfi^ A (x)4 B (x) + H.c. 

t "T^ V 7k 



(B4) 

Inserting Eqs. (fBl) in Eq. (fB4|) and Eq. l[2"Tjl. one even- 
tually obtains the estimate Eq. (f22|) for the impurity 
strength required to develop a resonance at w = in 
the site-like case. 

A similar expansion holds in the bond-like and in the 
hollow-like cases, Eqs. (|2"4"|) and ([23]) . respectively involv- 
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ing the constants 

A b = --^j2{ el{k ' S3 ~ 6k) + Rx ) ~°- 67 > ( B5a ) 

k 

Ah = jf 12 ( e2l9k7k + H - c -) ~ 2 - 35 ' (B5b ) 

k 

s * = iEw a = 3 > ( B5c ) 

k 

where e l9k is defined by Eq. (fTT|> . 

APPENDIX C: CURRENT DENSITY VECTOR 
WITHIN THE TIGHT-BINDING 
APPROXIMATION 

Here, we summarize some of the results employed to 
derive the expression of the conductivity in Sec. IIVBI 
within the tight binding approximation outlined in Ap- 
pendix [B] We begin by reminding the explicit expression 
of the Fourier transform of the paramagnetic component 
of the density current vector in reciprocal spaced 

JV ^ = -2^/^ 2q+k H^ (CD 

In the homogeneous limit (k = 0), one has^S 

J v (Q) = ^[H,r], (C2) 
in 

where H is the system's Hamiltonian, including the im- 
purity contribution. Exploiting Eq. (|C2D . one finds 

(ka|J v (0)|k'/3) = ie|«5 k ,, k+G e <G ^h a/3 (k), (C3) 

where a, (3 € {A, B}, G is a vector of the reciprocal lat- 
tice, and 8p = if j3 = A, and 8p = S3 if f3 = B. Due to 
the discrete translational invariance and the hcrmiticity 
of the current density operator, the adimcnsional matrix 
elements h a ^(k) fulfill the additional properties 

h AA {k) = h BB (k), (C4a) 
h AA (k) = -h AA (-k), (C4b) 
h AB (k) = -h BA (-k). (C4c) 



Moreover, the off-diagonal elements afford the explicit 
expression 



h AB (k) = V k 7k 

a 

1 3 

= -12 5 ^' S \ (C5) 

t=l 

to leading order in the overlap parameter s, where use 
has been made of Eq. |5]), which, together with Eq. (|C4cl) . 
yields the off-diagonal terms of the matrix elements. 

In order to find the diagonal terms, it is useful to ob- 
serve that the pseudoatomic wavefunctions introduced in 
AppendixlBlarc cylindrically symmetric. This implies the 
following overlap and dipole element for pseudoatomic 
wavefunctions centered on nearest neighbor sites 

J dr<l>(r)<f>(r±6t) = s, (C6a) 
J dr<Xr)r<Xr±<y = t\sS,, (C6b) 

where s is the band asymmetry parameter (Appendix [B| 
and I = 1, 2, 3, or, more compactly, 

J dr0(r)r0(r-r') = ^r' J dr <p(r) (f)(r - r'), (C6c) 

Making use of Eq. l|C6b[) . one eventually finds 



h QQ (k) = ---V k/ 5k 
a Z 

= Yal2 12 (Se-SmV k < S '- S -\(C7) 

i=l m=l,m^e 

with /3 k given by Eq. (|B2|) . 
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